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Abstract. We say that two graphs arc similar if their adjacency matrices are similar matrices. 

We show that the square grid Gn of order n is similar to the disjoint union of two copies of the 

(2) 

quartered Aztec diamond QADn—\ of order n — 1 with the path P„ on n vertices having edge 
weights equal to 2. Our proof is based on an explicit change of basis in the vector space on which the 
adjacency matrix acts. The arguments verifying that this change of basis works are combinatorial. 
It follows in particular that the characteristic polynomials of the above graphs satisfy the equality 
P(Gn) = P(Pi^') [P(QAZ)„_i)]^. On the one hand, this provides a combinatorial explanation for the 
"squarishness" of the characteristic polynomial of the square grid — i.e., that it is a perfect square, up 
to a factor of relatively small degree. On the other hand, as formulas for the characteristic polynomials 
of the path and the square grid are well known, our equality determines the characteristic polynomial 
of the quartered Aztec diamond. In turn, the latter allows computing the number of spanning trees 
of quartered Aztec diamonds. 

We present and analyze three more families of graphs that share the above described "linear 
squarishness" property of square grids: odd Aztec diamonds, mixed Aztec diamonds, and Aztec 
pillowcases — graphs obtained from two copies of an Aztec diamond by identifying the corresponding 
vertices on their convex hulls. 

We apply the above results to enumerate all the symmetry classes of spanning trees of the even 
Aztec diamonds, and all the symmetry classes not involving rotations of the spanning trees of odd 
and mixed Aztec diamonds. We also enumerate all but the base case of the symmetry classes of 
perfect matchings of odd square grids with the central vertex removed. In addition, we obtain a 
product formula for the number of spanning trees of Aztec pillowcases. 



Introduction 

The number of spanning trees of the Aztec diamond graph ADn (see Figure 2.1 for an illustra- 
tion of AD^) was shown by Knuth [6] to be given by a simple explicit product involving cosines. 
Stanley [9] posed then the problem of determining the number of spanning trees of the quartered 
Aztec diamonds (Figure 2.2 shows the quartered Aztec diamond of order five). The author found 
a solution in the fall of 1996, when he was a Postdoctoral Fellow at the Mathematical Sciences 
Research Institute in Berkeley (see also [3]). This solution is presented in Section 2, and was 
the starting point of the current paper. In the meanwhile, a different solution has been found by 
Richard Kenyon, James Propp and David Wilson [5, §6.8]. 

A useful observation in Cvetcovic et al. [4] (reproduced here as Theorem 2.2) allows one to 
deduce the number of spanning trees of certain planar graphs — the graphs considered in this 
paper included — from the characteristic polynomial of a slight modification of their dual. Our 
approach is to determine these characteristic polynomials by showing that the involved graphs 
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"reduce" to disjoint unions of graphs whose characteristic polynornials arc known. More precisely, 
given a graph whose characteristic polynomial we need to find, we provide a block diagonal matrix 
similar to its adjacency matrix, so that the diagonal blocks are adjacency matrices of graphs with 
known characteristic polynomials. Our solution is combinatorial, in that it accomplishes this by 
providing simple, conceptual changes of bases. 

We apply this approach four times. First, we decompose this way the square grids, thus 
obtaining the characteristic polynomial of quartered Aztec diamonds. Next, we decompose the 
odd and mixed Aztec diamonds, thereby obtaining the characteristic polynomials of the odd 
and mixed "halved" Aztec diamonds (see Section 3). In Section 6 we handle similarly the 
"Aztec pillowcase" graphs (we note that these are not related to the "Aztec pillows" previously 
introduced in the literature by James Propp). 

The results of Sections 2 and 3 allow us to solve most cases of two other natural enumeration 
questions: the symmetry classes of spanning trees of Azted diamonds (see Section 4), and the 
symmetry classes of perfect matchings of odd by odd square grids with the center vertex removed. 
We conclude the paper with a section posing some open problems. 

1. A similarity lemma for graphs possessing an automorphism of order two 
whose fixed points form a cut set 

The adjacency matrix of a weighted directed graph is the matrix whose rows and columns 
are indexed by the vertices and whose {u, t')-entry equals the weight of the directed edge from 
u to t; if there is such an edge, and otherwise. If the graph is undirected, one can replace 
each edge e by a pair of anti-parallel directed edges of the same weight as the weight of e, and 
use the previous definition. The characteristic polynomial of the graph G is the characteristic 
polynomial of its adjacency matrix; we denote it by P{G;x) (or simply P{G), if we need not 
display its argument). We think of unweighted graphs as weighted graphs in which all weights 
are equal to 1. 

Let G = {V, E) be a connected, undirected, weighted graph. For any subset X <zV, denote 
by {X) the weighted subgraph of G induced by X . 

Let T be an automorphism of G so that is the identity. Denote by Vq the set of vertices 
of G fixed by T. Assume that (y \ Vq) is the union of two disjoint graphs of the form (Vi) and 
(T(Fi)), for some suitable Vi CV (this happens whenever Vq is a cut set). 

Set G~^ := (Vi). Let G" :— {Vi U Vo)- Consider the directed graph on the vertex set Vi U Vq 
that has a directed edge from u to ii if and only if u and v are adjacent in G' . Weight this edge 
by twice the weight of the edge {u, v} if u G Vq and v G Vi, and by the weight of {u, v} in all 
remaining instances. Denote the resulting weighted directed graph by G~ . 

We say that the graphs Gi and G2 are similar — and write Gi ~ G2 — if their adjacency 
matrices are similar. 

Lemma 1.1. Under the above assumptions, G is similar to the disjoint union of G~^ with G~ . 

Proof. Let A be the adjacency matrix of G. For each vertex w of G consider an indeterminate By. 
Denote by U the complex vector space of formal linear combinations X^ugy c„e^, Cy € C. 

Let Fg be the linear map from U to itself that sends to ^^^^y ay^wCw, where Uy^w is the 
weight of the edge {v,w} if v and w are adjacent, and otherwise. Then the matrix of Fq in 
the basis B := {cy : v S V} is just the adjacency matrix of G. 

Set B' := {(e^ - eT(y))/2 :veVi} and B" := {{cy + eT(y))/2 -.veViU Vq}. One readily sees 
that B' U B" is a basis of U. 

2 



For any v gVi we have 



i^G((e.-eT(.))/2) = 



av,wew — (JT(j;),t«e 

= 2 — eT(io)) 

wev 

= 2 ^ «i),u;(ew — eT(«,))) 

where at the last equaUty we used V = ViUVoU T{Vi), the assumption that there are no edges 
between Vi and T{Vi), and the fact that the summand in the next to last line above vanishes 
when w e Vq. Thus B' spans an i^c-invariant subspace, and the matrix of its restriction to it is, 
in basis B', just the adjacency matrix of G~^. 
Similarly, for v G Vi we have 



i^G((e.+eT(.))/2) 



+ eT{w)) 



wev 



weVi weVo 

while for v €Vq 

Fciiev + eT(„))/2) = i?G(e^) = ^ a^,^e^ 

— ^ ^ ^VjW^w ~t~ ^ ^ ^v,w^w ^~ ^ ^ (^VjW^w 

weVi u)eT(Vi) tuevb 

weV\ weV\ weVo 

weV\ weVo 

Therefore, B" also spans an Fc-invariant subspace, and the matrix with respect to B" of the 
restriction of Fg to this Fc-invariant subspace is precisely the adjacency matrix of the weighted 
directed graph G~ . 

Thus the matrix of Fg in the basis B' U B" is block-diagonal with the two blocks equal to the 
adjacency matrices of G+ and G~ , respectively, and the claim follows. □ 

For the remainder of this paper, given a graph G = (V, E) we will denote by Ug the complex 
vector space of formal linear combinations X^^gyC^e^, c„ e C, the e^'s being indeterminates, 

3 



Figure 2.1. The Aztec diamond AD5. Figure 2.2. The quartered Aztec diamond QAD5. 
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Figure 2.3. The square grid G5. Figure 2.4. The weighted directed graph QAD^. 

and by Fq the linear map from Ug to itself that sends to J2wev where a„,^ is the 

{v, ti;)-entry of the adjacency matrix of G. 

2. The quartered Aztec diamond 

The Aztec diamond graph of order n, denoted ADn, is the subgraph of the grid (Z + 1/2)^ 

induced by the vertices {x,y) with \x\ + \y\ < n (Figure 2.1 shows AD^). The quartered Aztec 
diamond QADn is the subgraph of ADn induced by the vertices in its southeastern quarter 
{QAD5 is pictured in Figure 2.2). 

Denote by Gn the nx n grid graph (illustrated in Figure 2.3 for n = 5), and by P„ the path 
on n vertices having edge weights equal to 2. 

(2) 

Theorem 2.1. G„ is similar to the disjoint union of two copies of QADn-i with Pn ■ 

Gn ~ QADn-1 U QADn-l U P^^h (2.1) 



The proof is presented after Corollary 2.3 below. 

Theorem 2.1 implies in particular that the characteristic polynomial of Gn equals the product 
of the characteristic polynomials of the graphs on the right hand side of (2.1). 

It is well known (see e.g. [8, problem 1.29]) that the eigenvalues of the path P„ (with edge 
weights equal to 1) are 

J2C0S— -,2cos— -,...,2cos— -h 2.2 
[ n+1 n+1 n + lj 

the characteristic polynomial of P„ readily follows from this. 

On the other hand, the square grid G„ is the so-called tensor sum of P„ with itself (see e.g. 

[7]), and thus by [7, Theorem I] its eigenvalues are given by all possible sums of two of the 
numbers in (2.2): 

2cos^^ 1- 2 cos , l<j,k<n. (2.3) 

n+1 n+1 - V ^ 
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It follows that 



P{QADn_i;x) = n (a;-2cos^f- -2cos-^ ) . (2.4) 
-'--'■V n + 1 n+1 J 

l<i<j<n ^ ^ 

The following result of [4] provides a useful connection between characteristic polynomials and 
the number of spanning trees of certain planar graphs. For completeness, we include the proof. 
For a graph G we denote by t(G) the number of its spanning trees. 

Theorem 2.2 [4]. Let G be a connected planar graph all of whose bounded faces have r sides. 
Let G' be the graph obtained from the planar dual of G by deleting the vertex corresponding to 
the infinite face. Then 

t(G)=P(G';r). 

Proof. By a well known fact, G has the same number of spanning trees as its planar dual G""-. 
By the Matrix Tree Theorem (see e.g., [1]), t(G^) is equal to the determinant of the matrix 
obtained from the negative Laplacian of G^ (i.e., the difference between the diagonal matrix 
of its vertex degrees and its adjacency matrix) by deleting the row and column indexed by the 
vertex of G"*- corresponding to the infinite face of G. The hypothesis implies that the latter 
matrix is precisely the evaluation of the characteristic polynomial of G' at r. □ 

Corollary 2.3. The number of spanning trees ofQADn equals 



n 



l<j<k<n-l 



4-2cos:^ 2 cos — . (2.5) 

n n J 



Proof. This is a direct consequence of Equation (2.4), Theorem 2.2, and the fact that the graph 
obtained from the planar dual of QADn by deleting the vertex corresponding to the infinite face 
is isomorphic to QADn-2. □ 

Proof of Theorem 2.1. Apply Lemma 1.1 to the square grid G„, choosing T to be the symmetry 
across one of its diagonals (see Figure 2.3). The resulting graph G"*" is isomorphic to QADn-i. 
The resulting weighted directed graph G~ is obtained from QADn by replacing all its edges by 
pairs of anti-parallel dircx'ted edges of weight 1 , and then changing to 2 the weights of all directed 
edges originating from the n vertices on the hypotenuse of the convex hull of QADn, denote the 
latter by QAD„ (Figure 2.4 shows QADc^; edges between non-marked vertices mean that there 
are directed edges of weight 1 between both corresponding ordered pairs of vertices; an edge 
between an unmarked vertex v and a marked vertex w signifies a directed edge from to w of 
weight 1, and a directed edge from tz; to u of weight 2). We obtain 



G„ ~ QAD„-i U QADn- (2.6) 

Let U = ^ QAD ^^'^ ^ = ^qTd ^'^ defined at the end of Section 1. 

As proved in Lemma 2.6 below, the union of the set B of n(n — l)/2 vectors in the statement 
of Lemma 2.4 with the set B' of n vectors in the statement of Lemma 2.5 forms a basis of U. By 
Equation (2.9), the subspace spanned by B is F-invariant. Lemma 2.5 proves that the subspace 
spanned by B' is F-invariant as well. Thus the matrix of the linear map F in the basis B \J B' 
is a block diagonal matrix with two blocks. By (2.9), the first of these is precisely the adjacency 
matrix of QADn-i. By (2.13), the second block equals the adjacency matrix of the path P„ ' . 

Thus 

QADn ^ QADn-i U Pi^l (2-7) 

This and (2.6) imply (2.1). □ 



-> c 



4 
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Figure 2.5(a). The twelve cases. Figure 2.5(b). Twelve representatives. 

In order to present the lemmas employed in the above proof, we coordinatize the vertices of 



QAD^ using matrix-style coordinates, assigning the top left vertex coordinates (1, 1). Let be 
the indeterminate corresponding to vertex {i,j)mU = U q^j . Thus {cij > <n+l} 

is a basis of U, and the matrix of F = in this basis is the adjacency matrix of QAD^. 

Lemma 2.4. For i,j>l,i + j< n, define 



{if any index is outside the range [l,n], we omit the corresponding term). 
Then for all i,j>l,i + j< n, we have 



(2.8) 



(2.9) 



where by convention undefined Vki 's on the right hand side {these occur when k and I fail to 
satisfy both k,l >1 and k + I < n) are omitted. 



Proof. By the definition of F, we have 



(2.10) 



for all i, j > 1, i + j < n, where we omit the terms on the right hand side corresponding to 
undefined ek,i's, i.e. to those index values that fail to satisfy both k,l > 1 and k + I < n + 1. 

Call the union of the neighborhoods of all neighbors of a vertex in a graph the second neigh- 
borhood of that vertex. Then by equations (2.8), (2.10), and the hnearity of F, both sides of (2.9) 
are linear combinations of efc;'s corresponding to vertices (fc, I) from the second neighborhood of 
vertex {i,j). 

It is apparent from the definition oiQAD^ that, up to refiection in its symmetry axis, there are 
only twelve translationally distinct second neighborhoods of its vertices; the induced partition of 
the set of vertices is shown in Figure 2.5(a) (pairs of vertices of QAD^^ that arc mirror images with 
respect to its symmetry axis are understood to be in the same class of the partition). Because of 
this, it is enough to check that (2.9) holds for one representative of each of the twelve classes — for 
instance those indicated in Figure 2.5(b). 

These twelve checkings are accomplished by Figure 2.6; the order of the panels corresponds 
to the numbering of the representative vertices shown in Figure 2.5(b). 
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Figure 2.6. Checking the twelve cases. 
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In each panel there are two copies of the second neighborhood of the vertex in question. 
The vertex is indicated by the black dot. The top copy performs the computation of the 
left hand side of (2.9), using (2.8), (2.10), and the linearity of F. Consider for instance the top 
half of the sixth panel. By (2.8) and the linearity of F, 

F{vij) = F{ei-ij) - F{eij-i) + F{ei+ij) - F{eij+i). 

Each term ±F{ek,i) on the right hand side above is depicted in the figure by placing the corre- 
sponding uncircled plus or minus sign next to the vertex {k,l). By (2.10), we have 

F{ei-i,j) - F{eij-\) + i^(ei+ij) - F(ej,j+i) =(ei_2,j + ei_i,j_i + Cij + ej_ij+i) 

— (ei-i,j-i + ei,j_2 + ej+ij-i + ejj) 
+(6»,j + Cj+i j_i + ei+2,j + Cj+ij+i) 

— (ei-i,j+i + + ej+ij+i + j+2). 

The sixteen ±efc,;'s on the right hand side are marked on the figure by placing circled plus or 
minus signs next to vertex {k,l). Then the coefficient of any ek,i in F{vij) is visually apparent 
from the top half of the sixth panel: Simply "add up" the circled signs, counting a circled plus 
as 1, and a circled minus as —1 (the coefficient is clearly zero unless vertex [k, I) is shown in the 
figure). 

To continue with our example, turn now to the bottom copy of the sixth panel in Figure 2.6. 
It performs the computation of the right hand side of (2.9), in the following sense. By (2.8), we 
have 

+ Vij-i + Vi+ij + Vij+i ={ei-2,j — Ci-ij-i + Bij — ej-ij+i) 

+ (ej-l,i-l — 6i j_2 + ^i+l.j-l — 

+ {ei-l,j + l — SiJ + P-i+lJ + l — eiJ+2)- 

Each ±ek,i 's above is depicted in the figure by placing a circled plus or minus sign next to vertex 
{k,l). This affords a visual computation of the result, by the same adding of circled signs as 
above. 

One checks by inspection that in the top half of the sixth panel in Figure 2.6, each vertex gets 
the same sum of circled signs as the corresponding vertex in the bottom half. This establishes 
Equation (2.9) for vertices of type 6. 

The remaining cases are checked similarly by the rest of Figure 2.6. (The calculations required 
in panels 8 and 9 are identical to those in panels 5 and 6, respectively; thus we did not repeat 
them.) The double circled plusses and minuses in the last three panels indicate contributions of 
±2efe; , and appear due to the fact that the directed edge from a marked to an unmarked vertex of 
QAD^ has weight 2 (marked vertices are indicated by black diamonds). Note that in the bottom 
panels there are no contributions coming from the diamond vertices, because by the convention 
stated after Equation (2.9), the undefined terms on the right hand side of the latter are omitted 
(and by the statement after Equation (2.8), Vij is not defined if is a diamond vertex). □ 

Lemma 2.5. Let 

._f2, h3>l,i + j<n . . 

and define 

Wk := ^ Cijeij, (2.12) 

n-k+2<i+j<n+l 
i—jE{n—k,n—k—2,n—k—4,...,k—n} 
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Figure 2.7{i). Figure 2.7(m). Figure 2.7{iii). 
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Figure 2.7{iv). Figure 2.7{v). Figure 2.7{vi). 




Figure 2.7{x). Figure 2.7{xi). Figure 2.7{xii). 

for k = 1, ... ,71 {for n = 9 these are represented in Figures 2.7{i)-{ix)). Then 

F{wk)=2wk-i+2wk+i (2.13) 

for all k = 1, . . . ,n, where by convention on the right hand side we omit terms of index outside 
the range [l,n]. 
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Figure 2.8(a). The weight function wtg. Figure 2.8(b). The weight function wtio- 




Figure 2.9. d. 

Proof. This follows again by the linearity of F and the fact that its action on is the weighted 

sum of the indeterminatcs corresponding to the neighbors of vertex (i,j) in QAD,^. Figure 2.7 
corresponds to n = 9, but all the features needed to check the general case are present in it. Panels 
{i)-{ix) show the vectors w\, . . . ,wg, respectively: the vertices that contribute indeterminatcs 
with non-zero coefficients are indicated by the value of that coefficient next to them. 

Panel {x) shows ^F(w4) — one dot next to a vertex {k,l) represents the contribution of one 
unit to the coefficient of eu ■ It is apparent that this vector is the same as the sum of the vectors 
in panels {Hi) and (w). This checks (2.13) for fc = 4. 

Similarly, panels [ii), (iv), and (xi) check (2.13) for k = 3, while the case k = n = 9 follows 
by panels (viii) and {xii). □ 

Lemma 2.6. The vectors (2.8) and (2.12) form a basis o/U. 

Proof. Since there are n{n — l)/2 vectors Vij in (2.8) and n vectors Wk in (2.12), it suffices 
to show that eij is in their span for all i, j > 1, i + j < n + 1. We show first that the vectors 
Uk := Y,i+j=k+i eij, k = 1, . . . n, are in this span. 

Indeed, for 1 < fc < Z < n define 



[k,l] := 



(2fc-l) + (2fc + l) + --- + (2/-l), 

{2k) + {2k + 2) + --- + {2l), 



if n is odd, 
if n is even. 



Let wt„ be the weight function on the vertices of QAD^ whose nonzero values are given by the 
patterns shown in Figure 2.8 — Figure 2.8(a) shows the pattern for odd values of n, and Figure 
2.8(b) the pattern for even n (vertices of nonzero weight are marked by dots). 

Define the closed half-strip Ck to be the region described, in the standard rectangular coor- 
dinate system with origin at the northwest corner of QAD^, by {{x,y) : — (fc — 1) < a; + y < 
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Figure 2.10. (a). The sum in (2.14) for n and k odd. (b). Same for n odd and k even. 



k — l,x — y>k — 1} (for n = 9, C4 is illustrated in Figure 2.9). Color the vertices of QAD^ in 
chessboard fashion. Denote the set of vertices of QAD^ by V. 
We claim that for = 1, . . . , n — 2 one has^ 



E 



^in{x)vx-kwn-i-k + {k+'2)wn+i-k = (2n+4)(efe,i+efe_i,2H hei,fe) 



xeCfe+iny 

X of same color as vertex (k,l) 

(2.14) 

Indeed, suppose that n is odd. It is then readily seen that the definition (2.8) of the vectors 
Vij and the definition of the weight wt„ imply that when expressing the sum on the left hand 
side of (2.14) in terms of the vectors e^ , the nonzero coefficients are given by the pattern shown 
in Figure 2.10 (the two pictures illustrate the cases of odd and even k). On the other hand, 
when expressing the combination of the last two terms on the left hand side of (2.14) as a 
linear combination of the e^ 's, the nonzero coefficients are given by the patterns indicated by 
Figure 2.11. It is apparent from these two figures that (2.14) holds. The case of even n is justified 
analogously. 

By (2.14), the vectors wi, . . . , u„-2 are seen to be in the span of the vectors (2.8) and (2.12). 
Clearly, w„_i = and u„ = wi. Thus all Mfe's, fc = 1, ...,n are in the span of the vectors (2.8) 
and (2.12), as claimed at the beginning of this proof. 

We now show by induction on i + j > 2 that eij is in the span S of the vectors (2.8) and 

. . . , Un}. 

If i + i = 2, we must have i = j = 1, and the claim follows as en = ui. Assume that Cij is 
in S for all i + j < k. Then the definition (2.8) of Vk-ij implies that ek-i,i+i — Rfc-i+i.i is in 
S, for i = 1, . . . , A; — 1. Thus Ck-is+i = fife.i + Si with Si G S, for i — 1, . . . ^k — 1. Therefore 

Uk = efe,i+efc_i,2H Vei^k = fcefc.i + s, with s e S, implying ek,i G S. Since ek-i,i+i = ek^i + s^, 

it follows that Ck^i^i+i € 5 for all i = 0, . . . , A; — 1, and the induction step is complete. This 
completes the proof. □ 

Remark 2.7. Theorem 2.1 shows considerably more than the fact that the characteristic poly- 
nomial of the graph on the left hand side of (2.1) is equal to the product of the characteristic 
polynomials of the graphs on the right hand side of (2.1): it shows that the corresponding adja- 
cency matrices are similar (i.e., they have the same Jordan form). The same remark applies to 
the results of Section 3. 



^For a vertex x of QAD^ we write Vx for Vij, where are the coordinates of x. 
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Figure 2.11. (a). The sum of the last two terms on the left hand side of (2.14) for n and k odd. 

(b). Same for n odd and k even. 

3. The halved Aztec diamonds 

One way to define the Aztec diamond ADn is to say that it is the graph whose vertices are 
the white unit squares of a (2n + 1) x (2n + 1) chessboard with black corners, two vertices being 
adjacent if they correspond to white unit squares that share a corner. The analogous graph 
on the black unit squares of this chessboard is called the odd Aztec diamond of order n and is 
denoted ODn (Figure 3.1 shows OD5). 

Similar considerations on a 2n x 2n chessboard lead to graphs on the white and black unit 
squares that are isomorphic. We denote them by MD„ ("mixed" diamonds; MD5 is pictured in 
Figure 3.2). 

Define HODn to be the subgraph of ODn induced by the black unit squares on or above a 
diagonal of the (2n + 1) x (2n + 1) chessboard. Similarly, let HMDn be the subgraph of MDn 
induced by the black vertices on or above the black diagonal of the 2n x 2n board {HOD5 and 
HMD5 are shown in Figures 3.3 and 3.4, respectively). 

Let Qn be the graph obtained from the path P„ by including a loop of weight 2 at each vertex. 
Let Q'^ be graph obtained similarly from P„, but by weighting each loop by —2. The graph i?„ is 
obtained from Qn by changing the weight of the loop at the last vertex to 1, and R'^ is obtained 
from Q'^ by changing the weight of the loop at the last vertex to —1. 



Theorem 3.1. We have 



MDn ~ HMDn-l U HMDn-l Rr,\J R' 



(3.1) 



Proof. Apply Lemma 1.1 to MDn, choosing the automorphism T to be the reflection across 

its symmetry axis. The resulting graph G"*" is isomorphic to HMDn-i- The resulting weighted 
directed graph G~ , which we denote HMDn, is obtained from HMDn by: (1) marking the 2n 
bottommost vertices, (2) replacing each edge between two marked or two unmarked vertices by 
two anti-parallel directed edges of weight 1 , and (3) replacing each edge between an unmarked 
vertex v and a marked vertex w by a directed edge {v, w) weighted 1 and a directed edge {w, v) 
weighted 2 (Figure 3.9 illustrates this when n = 5). We obtain 



MDn ~ HMDn-l U HMDn- 



(3.2) 
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Figure 3.1. OD5. 



Figure 3.2. MD5. 



Figure 3.3. HOD5. 



Figure 3.4. HMD5. 
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Figure 3.5. Q5. 
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Figure 3.7. R5. 
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Figure 3.6. Q'^. 
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Figure 3.8. R'^. 




Figure 3.9. HMD5. 



Let t; be a vertex of HMDn- Proceed from v along a ray in the northeast direction. If no other 
vertex of HMDn is on this ray, set NE(i;) := 00; otherwise, let w be the first encountered vertex 
and set NE(t') := w. Define NW(u), SW(w), and SE(v) analogously, via the vertices closest to v 
in the northwestern, southwestern, and southeastern directions, respectively. 

For each vertex v of HMDn, let e„ be an indeterminate. Denote by Vi the set of unmarked 
vertices of HMDn- Let Bi C Vi consist of the n — 1 vertices on the northeastern side of the 
convex hull of Vi , and B2 C Vi of the n — 1 vertices on the northwestern side of the latter. Define 

' -eNE(i)) + eNW(i;) - ^SW{v) + eSE(i;)5 V G Vi \ {Bi U B2) 

fv ■= { ev + cnwCu) - esw(?;) + esE(?;), V G Bi (3.3) 



-Bv — eNE(«) - esw('u) + esE(v), 
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V G B2, 



Figure 3.10. The sets Sk for n = 6. 



where by definition Cqo := 0. It will turn out that these particular special definitions for the 
vertices in Bi and B2 are just the right ones for obtaining the desired decomposition of HMDn- 
Regard (3.3) as vectors in U = HMD r^' ^ ^^^^ analysis analogous to that in the proof of 
Lemma 2.4 shows that the action of F = Fjjj^^ on the vectors (3.3) is given by 

F{,fv) = /N(t)) + fw{v) + fsiv) + /e(u). (3-4) 

where N(i;), W(v), S{v), and E(f) are defined in analogy to NE(t)) — via rays from v in the 
directions of the four cardinal points — and foo '■= 0. 

The set V of vertices of HMDn naturally splits into n levels denote them, from bottom to 
top, hy Li, ... ,Ln- The same set can also be regarded as consisting of 2n columns of vertices — 
denote them, from left to right, by Ci, . . . , C2n- We define subsets Si, . . . ,Sn CV as follows. 

Color the vertices of V black and white in chessboard fashion so that Ci is black; let Vb be 
the set of black vertices. For 1 < k < n, act 

Sk := (Vb n (Ci U C2 U • • • U Cfe U C2n-k+i U C2n-k+2 U • • • U C2„)) 

U ((Cfc U Ck+i U • • • U C2n-k+i) n {Lk U Lk-2 U Lk-4 U • • • )) (3.5) 

(for n = 6, these arc pictured in Figure 3.10). 

Let wt be the weight function on V that assigns 1 to the vertices in Li, and (— l)''"^ • 2 to 

the vertices in Lk, if fc = 2, . . . ,n. Let wt'(u) := wt{v) if v is black, and wt'{v) := — wt(w) if v is 
white. 

Define, for A; = 1, . . . , n, 

Qk := (-!)*=-' J2 (3.6) 
veSk 

and 

9'k ■= wt'{v)ey. (3.7) 
veSk 

An argument similar to the one that proved Lemma 2.5 shows that the action of F on these 
vectors is given by 

F(.gi) = 2,gi+.92, (3.8) 

F{gk) = 9k^i + '2gk + gk+u 2<k<n-l, (3.9) 

F{gn) = gn-i + gn, (3.10) 



and 

Fig[) = -2g[+g'2, (3.11) 
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Figure 3.11. The support of an,k (circled vertices) and (dotted vertices). 



Fig'k) =9^-1-^ + 9^1^ 2<k<n-l, 



(3.12) 



F{9'n)=9n-l-9'n- 



(3.13) 



Furthermore, by Lemma 3.2, the span of the vectors (3.3), (3.6) and (3.7) contains Uk = 
'^veLk ^^'^ ~ ^veLk ^v^vi k = 1, . . . , n, where Cy equals 1 or —1 according as v is black 
or white. Then a simple inductive argument shows that each ey, v G V is contained in the span 
of the union of the vectors (3.3) with the u^'s and the w^'s. This implies that the union of the 
vectors (3.3), (3.6) and (3.7) forms a basis of U. 

However, by (3.4) and (3.8)-(3.13), the matrix of F in this basis is a block diagonal matrix 
consisting of three blocks. One of them, corresponding to the rows and columns indexed by the 
vectors (3.3), is by (3.4) the same as the adjacency matrix of HMDn-i- The other two, by (3.8)- 
(3.13), are the same as the adjacency matrices of i?„ and R'^, respectively. This implies (3.1). □ 



Lemma 3.2. The span of the vectors {fy : v G Vi}, {gk : k = 1, . . . , n}, and {g'j^ : k = 1, . . . ,n} 
given by (3.3), (3.6), and (3.7) contains Uk = '^y^^k^^ ^'^^ "fe ~ SweL^ c^^fj ^ = l,...,n, 
where Cy equals 1 or —1 according as v is black or white. 



Proof. Suppose n is even. For 2 < A; < n, let an,k be the function on V that is zero on 
the vertices not circled in Figure 3.11, and whose value at a circled vertex is given by the 
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corresponding entry in the following array (for briefness we denote a = n — k): 

k-l 

-(fc-1) -2(fe-l) 
fe-1 2(fe-l) 3(fe-l) 



fe-1 2(fe-l) 3(fe-l) ■•■ (a-l)(fe-l) 

-(fe-1) -2(fe-l) -3(fe-l) ••• -(o-l)(fe-l) -a(fe-l) 

fe-1 2(fc-l) 3(fc-l) ••• (a-l)(fe-l) a(fc-l) (a+l)(fe-l) 

fe-2 2(fe-2) 3(fc-2) ••• (a-l)(fe-2) a(fe-2) (a+l)(fe-2) 

1-2 2-2 3-2 ■■■ (a-l)-2 a-2 (a+l)-2 

12 3 ■■■ a-1 a o+l 

(read the pattern above from bottom to top; for k even the pattern ends on top as indicated; for 
k odd the alternation of signs along the rows of the top triangular portion of the array causes 
the entries in the top three rows above to be the negatives of the shown ones) . 

Let Pn,k be the function on V that is zero on the vertices not dotted in Figure 3.11, and 
defined at each dotted vertex v by I3n,k{v) = —an,k{v'), where v' is the reflection of v across the 
vertical symmetry axis of HMDn (note that the supports of an,k and (3n,k are disjoint). 

It is straightforward to check that for fc = 2, . . . , n one has 

X] i<^n,k{'") + l^n,k{v)) fv + [-{n -k + l)gk-2 - 9k-i + {n-k + 2)gk] = (2n + 1) X] 



Furthermore, if is plus or minus one according as v is black or white, one similarly checks that 

^ Cv (a„,fc(i;) + PnAv)) fv + {-if [{n-k + l)g'k_2 + g'k-i -{n-k + 2)gi] 
veVi 

= (2n + 1) ^ c^Gy, 
veLk 

for = 2, . . . , n. Since by definition ui = gi and u[ = g[, the above equalities prove the statement 
for n even. 

The above equalities hold without change also for odd n, provided we define an,k and Pn,k to 
be the negatives of their values above. □ 
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Figure 3.12. HODi^. 



Theorem 3.3. We have 

ODn ~ HODn-i U HODn-1 U Pi U g„ U g;. (3.14) 

Proof. We proceed in a way analogous to the proof of Theorem 3.1. Lemma 1.1 imphes 



0£)„ ~ HODr,-i U HODr^, (3.15) 



where HODn is the weighted directed graph obtained from HODn by marking its bottommost 
2n + 1 vertices, replacing each edge of it between two unmarked or two marked vertices by a pair 
of opposite arcs of weights 1, and each edge between an unmarked vertex v and a marked vertex 
w by an arc {v,w) of weight 1 and an arc {w,v) of weight 2 (Figure 3.12 shows HOD4). 

Let V be the set of vertices of HODn, and let V\ be the set of its unmarked vertices. For each 

V V, let e-u be an indeterminate. For v Vi, define f-^ by (3.3), the same formulas we used 
in the proof of Theorem 3.1 — except for the case when v is the topmost vertex in Vi, when we 
define := — esw(t>) + '^se(«)- Then, as an argument similar to the one in the proof of Theorem 
3.1 readily checks, (3.4) holds for all v E Vi. 

The subsets of vertices and the weights on V that will provide us with the vectors we need 
to complete {/t, : t; e Vi} to a basis of U = ^ hod „ ^'^^ defined slightly differently than 
in the proof of Theorem 3.1. Define the g^'s and g^'s to be the linear combinations of the e„'s 
indicated by the patterns shown in Figure 3.13(a) and (b), respectively (these figures correspond 
to n = 4). 

An argument similar to the one given in Lemma 2.5 shows that the action of F = Fjjq^ on 
these vectors is given by 

F{gr)=2gi+g2, (3.16) 
F{gk) = gk-i + 2gk + gk+i, 2<fc<n-l, (3.17) 
F{gn)=gn-i+2gn, (3.18) 

and 

F{g[) = -2g[+g'^, (3.19) 

F{9'k)=9'k-i-^9k+9'k+i^ 2<k<n-l, (3.20) 

F{9'J^g'n^,-2g'n. (3.21) 

Finally, the vector h that equals the alternating sum of e„'s over the set of the n vertices of 
HODn on the northwestern side of the convex hull of its vertex set (the topmost e„ having 
coefficient +1) is readily seen to satisfy F{h) = 0. 

An inductive argument similar to the one used in the proof of Lemma 2.6 shows that for all 

V &V, By is in the span of the vectors {/„ : v gVi}, {uk : fc = 1, . . . , n}, {u^ : k = 1, . . . ,n} and h 
(the only part that requires separate justification is the checking of the base case; this is provided 
by Lemma 3.5). Thus, by Lemma 3.4, the span of the vectors {/„ : v G Vi}, {gk : k = 1, . . . , n}, 
{g'k : k = 1, . . . ,n}, and h contains e„ for all v €V, and hence these vectors form a basis of U. 
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Figure 3.13(a). The ^^fc's for n = 4. 
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Figure 3.13(b). The g'^'s for n = 4. 



By the above formulas describing the action of F on them it follows that the matrix of F in 
this basis is a block-diagonal matrix consisting of four blocks: one the same as the adjacency 
matrix of HODn-i, by the fact that (3.4) holds for V\] the next two equal to the adjacency 
matrices of the path-like graphs Qn and Q^, respectively, by (3.16)-(3.21); and the last equal to 
the 1x1 block consisting of a single 0, i.e., the adjacency matrix of the path Pi. Together with 
(3.15) this implies (3.14). □ 



Lemma 3.4. The span of the vectors {/„ : v £ Vi}, {gk : k = 1, . . . , n}, {g^ : fc = 1, . . . , n}, and 
h in the proof of Theorem 3.3 contains the vectors Uk and u'j, defined by the patterns in Figure 
3.14, fork=l,...,n. 



Proof. Let wt„ be the function on the vertices of of HODn given by the pattern 
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Figure 3.14(a). The UkS for n = 4. 
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Figure 3.14(b). The u'^'s for n = 4. 





-(n-l) n-1 
-2(n-2) -(n-2) n-2 2(ti-2) 

-3{n-3) ■■■ -6 -3 3 6 ••• 3(n-3) 

-2(n-2) ••• -6 -4 -2 2 4 6 ••• 2(n-2) 

-(n-1) ■■■ -4 -3 -2 -1 1 2 3 4 ••• n-1 

O--- 000 000000 0---0 

Let Cfe be the closed shaded region illustrated in Figure 3.15. It is straightforward to check 
that 

^ fftn{v)fy + [-{n -k + l)gk-i + {n-k + 2)gk] = (n + l)^^, 
veViDCk 
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Figure 3.15. The region d for HO Dr. 

for all = 1, . . . , n. 

Similarly, if wt', is the weight funetion on the vertices of HODn that equals wt„ on the black 
vertices and — wt„ on the white vertices, one checks that 

^ wt'„{v)U + (-1)^-1 [{n-k + + {n-k + 2)gi] ={n+ 1)4, 

veVinCk 

for A; = 1, . . . , n. □ 

Lemma 3.5. Let t be the topm,ost vertex of HODn- Then Ct is in the span of the vectors 
{fv ■ V G Vi}, {iik : fc = 1, . . . ,n}, {ttj. : fc = 1, . . . ,n}, and h. 



Proof. Let a„ be the weight function on the subset Vi of the vertices of HODn obtained by 
weighting the ith topmost row, from left to right, by the successive coordinates of the vector 

(-l)^(n ^ + 1) _ 1, 0, 2i - 3, 0, . . . , 3, 0, 1), 

2t 

for i = 1, . . . , n (recall that by definition HODn has n + 1 rows). 
Let 

(_l)"+fe(2A;(n - fc + 2) - n - 1) 
4(n-fc + l)(n-fc + 2) ' 

Then it is straightforward to check that 

n 

X] <^n{v)fv + ^dk (ttfc + u'k) + h= {n + l)et. 
veVi fe=i 

This proves the claim. □ 
Corollary 3.6. The number of spanning trees of the halved Aztec diamond HMDn is given by 

t{HMDn)= n P~^^°®2^^°^S) ■ ^^"^^^ 

l<j<k<2n-l ^ ^ 
j+k<2n-l 

Proof. The graph obtained from the planar dual of HMDn by removing the vertex corre- 
sponding to the infinite face is readily seen to be isomorphic to HODn-2 (see Figures 3.3 and 
3.4). Thus, Theorem 2.2 implies 



t{HMDn) = P{HODn-2-A). 
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(3.23) 



By [6] we have 

P(02.„;.)=.n(.-4cos^^) n (-^cos^cos^). (3.24) 

i=l ^ ^ l<i<fe<2n+l ^ ^ 

By (2.2), the eigenvalues of Q„ are 2 + 2cos^^ = 4cos^ 2ii+2 • ^ ~ 1, . . . ,n, while those of 
Q'^ are —2 + 2 cos j — 1, . . . , n; the latter hsted in reversed order thus equal —4 cos^ 2nT2 ' 
J = 1, . . . , n. Theorem 3.3 and (3.24) imply then, after some manipulation, that 

l<j<k<2n+l ^ ^ 
j+fe<2n+l 

The statement follows now by (3.23). □ 
Corollary 3.7. The number of spanning trees of the halved Aztec diamond HODn is given by 

t{HOD^)= n (4-4-^2;^ -^2^)- (3.26) 

l<j<k<2n ^ ^ 
j+k<2n 

Proof. The graph obtained from the planar dual of HODn by removing the vertex correspond- 
ing to the infinite face is readily seen to be isomorphic to HMDn-i- Therefore, by Theorem 2.2 
we obtain 

t{HODn) = P{HMDn-i-A). (3.27) 

By [6] we have 

P(MZ.„;.) = n(.-4cos^^) n (--^-2;7^-2;^^)- (^.28) 

i=l ^ ^ l<j<k<2n ^ ' 

It is not hard to show that the eigenvalues of the adjacency matrix of i?„ are 2 + 2 cos 2n+\ ~ 
4cos^ 2ra+i ' ^ ~ l)---)^) while the eigenvalues of the adjacency matrix of R!^ are just the 
negatives of the former. Using this, one sees after some manipulation that 

P{HMDr,-i;x)= J] ('x-4cos^^cos^-^y (3.29) 

j+k<2n 

Together with (3.27), this implies (3.26). □ 

Remark 3.8. It is amusing to note that the range for {j, k) in the product giving the character- 
istic polynomial of HMDn can naturally be regarded as the vertex set of HODn, and vice- versa. 

4. Symmetry classes of spanning trees of Aztec diamonds 

It is easy to see that no spanning tree of an Aztec diamond AD^ can be symmetric with 
respect to a a symmetry axis i of the diamond that makes a 45° angle with the horizontal. 
Indeed, suppose ADn had a spanning tree T symmetric about £. Let a and b be two distinct 
vertices of ADn on £. Since T is connected, there exists a path P in T connecting a to b. Since 
T is symmetric about £, the reflection P' of P across i is also contained in T. But no path in 
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Figure 4.1. A^. 



Figure 4.2. B5. 



Figure 4.3. Cs- Figure 4.4. Dg. 

ADn from a to & is invariant under this reflection, so P and P' are two distinct paths connecting 
two vertices of the tree T, a contradiction. 

The same argument works for the odd Aztec diamonds ODn when n > 2. The mixed Aztec 
diamonds MDn do not possess a symmetry axis at a 45° angle from the horizontal. 

Let h, V and r denote the symmetry across the horizontal, the symmetry across the vertical, 
and the rotation by 90°, respectively; view them as elements of the symmetry group of ADn 
or ODn- The above two paragraphs imply that there are a total of five inequivalent symmetry 
classes of spanning trees of ADn and ODn'. the base case — trees with no symmetry require- 
ment; horizontally symmetric trees — invariant under the action of h; horizontally and vertically 
symmetric trees invariant under the action of {h,v); trees invariant under the action of r^; and 
trees invariant under the action of r. Since MDn is not r-invariant, there are only four symmetry 
classes for its spanning trees. 

We denote the number of spanning trees of a graph G that are invariant under the action of 
the group of symmetries H by t_ff (G). 

The base case was done by Knuth [6]. In this section we provide product formulas for all 
but three of the remaining cases. The latter are phrased as open problems. Our arguments will 
express these numbers in terms of the number of perfect matchings of three families of subgraphs 
of the infinite grid 1? . Explicit formulas for the latter are deduced using results from the previous 
sections and the factorization theorem for perfect matchings of [2] . 

To define our families of graphs, it will be convenient to denote by NE{i,j) the infinite zig-zag 
lattice path in starting at the lattice point and taking alternately two unit steps north 
and two east; NW{i,j) is defined analogously, alternating between two steps north and two steps 
west. 

Let G2n be the subgraph of the infinite grid induced by the vertices with < i,j < 

2n — 1. We define An to be the subgraph of G2n induced by its vertices on or above the path 
NE{1,0). Bn is defined analogously, but using the path NE{0,0) instead of NE{0,1) {A5 and 
B5 are illustrated in Figures 4.1 and 4.2, respectively; the two dots in the latter emphasize that 
there are two vertices at the indicated positions). 

Define C„ to be the subgraph of induced by the lattice points (i, j), j > 0, that are on or 
below both NE{-n + 1, -1) and NW{n - 1, -1); Cg is shown in Figure 4.3. 
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A fourth family will be relevant in the proof of (4.3); it will also bo used in Section 5. We 
define £>„ to be the subgraph of 1? induced by the lattice points (i, j), j > 0, that are on or 
below both NE{-n, -2) and NW{n - 1, -1); is shown in Figure 4.4. 

For a graph G we denote the number of its perfect matchings by M(G) (if G is weighted, M(G') 
denotes the sum of the weights of all its perfect matchings, the weight of a matching being the 
product of the weights of its constituent edges). 



Lemma 4.1. We have 



and 



M{An) = ^ TT f4-2cos^^ -2cos-^V (4.1) 

l<j<k<n 



l<j<k<n ^ ' 

M(C„) = -^ n (4-2cos^-2cos^), (4.3) 

j-\-k even 

M(I.„) = ^ n (4-2cos^-2cos^). (4.4) 

j+k odd 



The following result of Temperley [10] will be useful in our proof. Let G = {V, E) be a finite 
connected planar graph; denote the set of its bounded faces by F. Pick a point v{e) in the 
interior of each edge e of G; let Ve := {v{e) : e G E}. Pick also a point v{f) in the interior of 
each bounded face / of G, and let Vf := {?'(/) : / G F}. Define T(G) to be the graph whose 
vertex set is VUVeUVf, with edge set consisting of the pairs {v, v{e)}, v GV, e G E, v incident 
to e, and {^^(e), v{f)}, e G E, f G F, e incident to /. In particular, if G„ is the nxn grid graph, 
T(G„) is the refined grid G2n-i. 

Theorem 4.2 (Temperley [10] [8, Problem 4.30]). For any v G V incident to the infinite 
face of G we have 

t(G) = M(T(G) \u). 



We will also employ the following "Factorization Theorem" we obtained in [2] . 

Theorem 4.3 (Factorization Theorem [2, Theorem 2.1]). Suppose the weighted bipartite 
planar graph G is drawn on the plane so that it is symmetric with respect to a horizontal symmetry 
axis £. Label the vertices on £ from left to right by ai,bi,a2,b2, . ■ . ,ak,bk, and assume they form 
a cut set. Properly color the vertices of G white and black, so that ai is white. 
Let G' be the graph obtained from G by 

(i) removing the edges above i incident to any white ai or any black bj; 

(a) removing the edges below £ incident to any black ai or any white bj; and 

{Hi) reducing by 1/2 the weight of any edge of G lying along £. 

Let G"*" be the subgraph of G' induced by the vertices above £, the black Oi 's, and white bj 's. 
Let G~ be the subgraph of G' induced by the vertices below £, the white a, '5, and black bj 's. Then 



M(G) = 2*=M(G+)M(G-). 
23 



Figure 4.5. 



Figure 4.6. 



Proof of Lemma 4-1- Apply the factorization theorem for perfect matchings [2, Theorem 2.1] 
to the graph H obtained from the square grid G2n+i by removing its top right vertex v, with 
respect to the symmetry axis of H (see Figure 4.5). One of the resulting subgraphs of H is then 
isomorphic to An (the boundary of this subgraph is indicated by the upper bold outline in Figure 
4.5), while the other can be regarded as T{QADn+i) \ v (the boundary of the latter is traced 
out by the lower bold outline in Figure 4.5). Since H itself can be regarded as T(G„+i) \ v, we 
obtain by the factorization theorem that 

M(T(G„+i) \v)=T M{An) M(T(Q^i?„+i) \ v). (4.5) 

By Theorem 4.2, the left hand side equals t(G„+i), which in turn by Theorem 2.2 equals 
P(G„;4); (2.3) provides an explicit formula for the latter. On the other hand, by Theorem 4.2, 
yi{T{QADn+i) \ v) = t(Q^D„+i), for which Equation (2.4) gives a product expression. Solving 
for M(A„) in (4.5) yields then (4.1). 

Equality (4.2) follows by applying the factorization theorem to the square grid G^n and its 
diagonal symmetry axis. The resulting subgraphs are both isomorphic to B„. We obtain 

M(G2„) = 2" [M(B„)]'- (4-6) 
A classical result of Kasteleyn, Temperley and Fisher (see e.g. [8, Problem 4.29]) gives 

l<j,fe<n ^ 

Plugging this into (4.6) and using the well known identity 11^=1 2 cos = 1 we obtain (4.2). 

Next, we prove (4.4). Apply the factorization theorem to the graph T(<5^-D„+i) \ v, with v 
chosen to be both on the symmetry axis and on the infinite face of T{QADn+i) (see Figure 4.6). 
One of the resulting subgraphs is then while the other is T(iJ0Z)„/2) \ when n is even, 
and T(iJM£)(„_|_i)/2) \ v, when n is odd. The formulas provided by Corollaries 2.3, 3.3 and 3.4 
can then be used to obtain an expression for M(_D„ ). After some manipulation one obtains (4.4). 

Now M(G„) can be worked out by applying the factorization theorem to the graph An- It is 
readily seen that one of the resulting subgraphs is while the other is G„. We obtain 

M(^„) = 2L^Jm(£)„)M(G„). 



Using (4.1) and (4.4) this implies (4.3). 
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Theorem 4.4. The only nonempty non-trivial symmetry class of spanning trees of ADn is that 
of horizontally symmetric trees. We have 



it^f,){ADn)=2ni{HMD^)=2n \{ U - 4 cos cos ^ j . (4.E 

l<7<fe<2n-l ^ ^ 



i<j<k<2n- 
j+k<2n-l 



Proof. Suppose T is a horizontally and vertically symmetric spanning connected subgraph of 
ADn. Since T is connected, it contains some edge e crossed by the horizontal symmetry axis 
Ih- The mirror image e' of e across the vertical symmetry axis £y is also in T, as T is vertically 
symmetric. Let a be the top vertex of e, and a' its mirror image across £„. Since T is connected, 
it contains a path P from a to a'. But then the union of the reflection of P across £h with the 
edges e and e' yields another path in T, distinct from P, connecting a to a'; so T cannot be tree. 

A similar argument shows that, if T is a spanning connected subgraph of AD„ which is 
invariant under rotation by 180°, and a is say one of the four centralmost vertices of AD^, then 
T must contain two distinct paths connecting a to its image through rotation by 180°. Since 
quarter turn invariant spanning trees are obviously also half turn invariant, to conclude the proof 
it suffices to prove (4.8). 

Let T be a horizontally symmetric spanning tree of ADn- The argument in the first paragraph 
of this proof implies that T contains precisely one edge e that crosses ih- The removal of this 
edge from T leaves two mirror images of a spanning tree of HMDn, and defines a 2n-to-l map 
from the set of horizontally symmetric spanning trees of ADn to the set of spanning trees of 
HMDn. This proves the first equality of (4.8); the second follows by Corollary 3.6. □ 

Theorem 4.5. The non-trivial symmetry classes of spanning trees of ODn not involving rota- 
tions are enumerated by 



^,){0Dn)=M{C2n-i) = ^ n (^-2^°'£-2cos^) (4-9) 



l<j<k<2n 
j+k even 



and 

tih,v){ODn) = M{An-i) = ^ H [4-2cos^-2cos— y (4.10) 



Proof. Let T be an /i-invariant spanning tree of ODn, and let a and b be the two vertices of 
ODn that are both on its infinite face and on its horizontal symmetry axis ih- Let P be the 
path in T connecting a to b. By the uniqueness of such a path it follows that P is /i-invariant. 
Clearly, this implies that P consists of the 2n edges of ODn along £h . 

It follows that the /i-invariant spanning trees of ODn are in bijection with the spanning trees 
of HODn that contain all the bottom 2n edges. In turn, the latter can be identified with the 
spanning trees of the graph G obtained from HODn-i by including a new vertex w connected 
to its 2n — 1 bottommost vertices. However, T(G) \'w \s readily seen to be isomorphic to C2n-\- 
Then (4.9) follows from Theorem 4.2 and (4.3). 

An argument similar to the one in the first paragraph of this proof shows that the (/i, v)- 
invariant spanning trees can be identified with the spanning trees of QADn which contain all the 
2n edges along the straight line portions of its boundary. In analogy to the previous paragraph, 
by Theorem 4.2 the number of the latter is seen to equal M(A„_i). Thus (4.10) follows by (4.1). □ 
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Theorem 4.6. The non-trivial symmetry classes of spanning trees of MDn not involving rota- 
tions are enumerated by 



twiMD,,) = M{C,n-2) = ^ n (4-2c°^2i^-2cos2^) (4.11) 



l<j<k<2n-2 
j+k even 



and 

n-l 



t(^,„)(MD„) = M(B„_i)=2"-incos-^ n (4cos2-^ + 4cos2 ^ 



2n- 1 



j=l l<j<fe<n-l 

(4.12) 

Proof. Arguments analogous to the ones in the proof of Theorem 4.5 show that t(/j)(M£'„) = 
M(C2„-2), and t(^h,v}{MDn) = M(S„_i). Then (4.11) and (4.12) follow by (4.3) and (4.2). □ 

Open Problem 4.7. Find formulas for t(^r){ODn), t(r2)(0-D„), and t(r2)(MZ)„). 

5. Symmetry classes of perfect matchings of odd squares with a central unit hole 

The existence of formula (4.7) invites one to look for an analog in the case of odd square grids. 
One very natural candidate is the graph obtained from the (2n + 1) x (2n + 1) square grid 
by removing its central vertex (see Figure 5.1). 

As it tiirns out, the numbers M(i?„) do not seem to have many factors in their prime factor- 
ization (see Section 7) . This makes it implausible for them to possess a "nice" product expression 
(e.g. similar in style to the product formulas that occur in this paper). 

However, all non-trivial symmetry classes of perfect matchings of Hn turn out to be enumerated 
by variations of the product formulas we have encountered in the previous sections. 

To present our results we will need, besides the graphs introduced in the previous section, to 
enumerate the perfect matchings of certain weighted versions of An and _B„. 

Define An to be the weighted graph obtained from An by weighting its top 2n — 1 edges by 
1/2, and keeping weight 1 for all its other edges. Define Bn to be the weighted graph obtained 
from Bn by weighting its top 2n — 2 edges by 1/2, and keeping weight 1 for all its other edges. 
For a weighted graph G, M{G) denotes the sum of the weights of all its perfect matchings, the 
weight of a matching being equal to the product of the weights on its constituent edges. 

Lemma 5.1. 

M(i„) = -LfT('4_4cos-^) n f4-2cos-^-2cos-^V (5.1) 

and 

nB.) = ^ n (4-4cos|^cos^). (5.2) 

j+k<2n-l, j-l-feodd 

Before giving the proof of the above lemma, we show how it implies the following result. 

Clearly, no perfect matching of Hn can be symmetric with respect to the diagonal. Thus the 
reflection h across the horizontal, the reflection v across the vertical, and the rotation r by 90° 
generate the relevant symmetry group. As is Section 4, there are four inequivalent non-trivial 
symmetry classes, corresponding to the subgroups (h), {h,v), (r^), and (r), respectively. Denote 
by Mk{G) the number of perfect matchings of G that are invariant under the symmetry group K. 
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Figure 5.1. Factorization applied to H^Kr"^). 



Theorem 5.2. We have 

ri — 1 2n 



M<.) = n n (4 - 2 cos ^|±i^ - 2 COS ^) , (5.3) 

j=0 fc=0 ^ ' 

M(.,>(H,„)= n (4^°^'^+4cos^^)> (5-4) 

(ir„)=2"-4tJ rr ('4-2cos^-2cos-^V> (5-5) 



'i-<j<k<n 
j-\-k odd 



M<.>(F2n-i)=2 n (4-2cos^-2cos^) [] " £ ' 



l<J<'s<"-l l<j<fe<2n-l 

j+/s<2n-l, j+feodd 



(5.6) 



M<.)(H2„) = ^n(4-4cos^) 

.7=1 ^ 



3-- 

X ^ 



n( 4 - 2 cos — 2 cos I ( 4 cos^ h 4 cos' 

\ '? n-\-^ ^/n -\-^ \ '? n-\-^ 



i<r<-k<n V 2n + 1 2n + 1; V 2n + 1 2n + 1, 

Proof. For a perfect matching of Hn to be /i-invariant, all vertices along the horizontal sym- 
metry axis ih must be matched by edges along £h; this is possible only when n is even. The 
/i-invariant perfect matchings of can then be identified with perfect matchings of the rectan- 
gular 2n X (4n-|- 1) grid graph. The latter were enumerated by Kasteleyn, Fisher and Temper ley; 
(5.3) follows from the form given in [5]. 

By a similar argument, the {h, w)-invariant perfect matchings of i?2n can be identified with 
perfect matchings of the square grid graph G2n- Thus (5.4) follows by (4.7). 

To prove (5.5), note that Mf^r^){Hn) = M(il„/(r^)), where Hn/{r^) is the orbit graph of the 
action of on i?„. This orbit graph can be regarded as being obtained from the subgraph of 
Hn induced by its vertices on or under a diagonal d by identifying the pairs of vertices on d that 
are at the same distance from the center; Figure 5.1 shows HQ/{r^). 
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Figure 5.2. Factorization for ifr/ (r) . 
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Figure 5.3. Factorization for ff6/(?')- 



Apply the factorization theorem of [2] to iJ„/(r^) with respect to the diagonal perpendicular 
to d. The resulting graphs both turn out to be isomorphic to 2?„. We obtain 

M<,2>(F„) = 2"[M(i)„)]2. 

Formula (4.4) implies then (5.5). 

Wc now turn to the proof of (5.6). Let _ff„/(r) be the orbit graph of the action of r on 
By definition, M(r)(-^^n) = M(i?„/(r)). The graph Hn/{r) is obtained from the subgraph of 
induced by its vertices on or under the union of its two diagonals by identifying the pairs of 
vertices on this union that arc at the same distance from the center; Figure 5.2 shows Hj/{r). 

Apply the factorization theorem of [2] to (r) with respect to the vertical symmetry axis. 
The shape of the resulting subgraphs is different for different parities of n. 

For odd n, one of the graphs is ^(n-i)/2, while the other is B(^n+i)/2 (see Figure 5.2). We 
obtain 

M(,)(iJ2n-i) = 22"-iM(A„_i)M(B„). 

Using (4.1) and Lemma 5.1, this implies (5.6). 

On the other hand, when n is even, one of the subgraphs obtained by applying the factorization 
theorem to Hn/{r) is -B„/2, while the other is An/2 (see Figure 5.3). Thus 

M(,)(/f2„)=22"M(B„)M(i„). 

Formula (5.7) follows then by (4.2) and Lemma 5.1. □ 

Proof of Lemma 5.1. Apply the factorization theorem of [2] to the graph C2n (see Figure 5.4; 
the dotted edges are weighted 1/2). One of the resulting subgraphs is then A„, while the other 
is Bn- We obtain 

M(C2„) = 2"M(i„)M(S„). 
Plug in the expressions for M(i?„) and M(C2n) given by (4.2) and (4.3). We obtain 

If 4 - 2 COS — 2 COS ■ 



22" 11 V 2n+l 2n+l 



j+k even 



2"M(i„) TT (4cos2-i^+4cos2 



\ 2n+l 2n+l/ /cox 

l<j<k<n ^ ' (5.8) 
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Figure 5.4. Factorization applied to Cs- Figure 5.5. Factorization applied to G5. 



Use the identity 2cos^(a;) = 1 + cos(2a;) to write 



4 cos^ — - + 4 cos^ = 4 + 2 cos ^ ^ + 2 cos ■ 



2n + l 2n+l 2n+l 2n+l 

, „ (2n + 1 - 27)7r „ (2n + 1 - 2k)n 

= 4 - 2 cos 2 cos '—. 

2n + 1 2n + 1 

Rewrite the factors on the right hand side of (5.8) using the above equality, and simplify out 
their occurrences in the product on the left hand side of (5.8). Solving for M(A„) in the resulting 

equation yields (5.1). 

To deduce (5.2), consider the graph G„ obtained from T{H M Dn) by including a new vertex 
V on its vertical symmetry axis i and joining v by an edge to the two bottommost vertices of 
HMDn on its infinite face (see Figure 5.5). Apply the factorization theorem of [2] to the graph 
G with respect to £. One of the resulting subgraphs is then while the other is T{QADn) \ w, 
where w is a leaf of QADn- We obtain 

M(G„) = 2" M(B„) M(T(gAI>„) \ w). (5.9) 

However, since v can be matched in G„ only to tu or its mirror image across we have by 
symmetry M(G„) = 2M(T(iJM£)„) \ w). Using Theorem 4.2, (5.9) becomes 

2t{HMDn) = 2^M{Bn)t{QADn). 
Combining this with formulas (2.5) and (3.22) we obtain 

n (4-4cosgcos^)=2"M(B„) n (4-2cos^-2cos^). (5.10) 

l<i<fe<2ri-l ^ ^ l<j<k<n-l ^ ' 

j+k<2n-l 

For same parity indices j and k in the product on the left hand side, write the corresponding 
factor as 

jw kn , „ (k — 7)7r „ (k + j)n 

4-4COS — cos — = 4-2cos-5^ ^ 2cos ^ . 

2n 2n 2n 2n 

The resulting factors cancel all the factors in the product on the right hand side of (5.10). Solving 
for M{Bn) in the resulting equality yields (5.2). 

□ 

29 



7 




4 


8 












3 


9 
















2 
























































12 
















19 


13 












18 


14 




17 



1 IC 
20 11 



7 




4 


8 












3 


9 
















2 
























































12 
















19 


13 












18 


14 




17 



1 

20 



15 16 15 16 

Figure 6.1. The Aztec pillowcase graph AP^. 



6. Aztec pillowcases 



Theorem 2.1 shows that the square grid graph G„ is similar to the disjoint union of two copies 
of the same graph (namely, QADn-i) with a graph having a relatively small number of vertices 
(namely, P^^^). 

Analogous statements about the Aztec diamond graphs MDn and OZ)„ follow by Theorems 3.1 
and 3.2. 

These three results are instances of finite graphs G with an automorphism T of order 2 so that 
if Gi is the subgraph of the orbit graph G/T induced by the orbits of the vertices not fixed by 
T, then G is similar to the disjoint union of two copies of Gi with some simple weighted graph 
on the vertices fixed by T. If in addition the number of vertices fixed by T is relatively small 
compared to the total number of vertices, we refer to such a pair {G, T) as linearly squarish. 

Which pairs {G,T) are linearly squarish? 

This seems to be a difficult question to answer. Lemma 1.1 applies whenever = 1 and the 
fixed points of T form a cut set. Its proof shows that G is similar to the disjoint union of G\ 
with a directed weighted graph constructed from the orbit graph G/T. The challenging part is 
to decide when the latter is similar to the disjoint union of Gi with some other graph. 

The question seems hard even when we restrict to the special case of graphs that are dual 
to simply connected lattice regions in Z^. If one uses the natural choice (2.8) (or (3.3)) to 
define part of the desired new basis, the required action (2.9) (or (3.4)) of F will hold as long as 
the corresponding vertex has its entire Z^-second neighborhood contained in G; but in general 
problems arise near the boundary. We may regard the above examples of the square grid and 
mixed and odd Aztec diamonds as instances when we managed to solve these problems at the 
boundary. 

In this section we give another example of a linearly squarish graph, whose being so is based 
on a variation of the change of basis (2.8). 

Define the Aztec pillowcase graph APn to be the graph obtained from two copies of ADn by 
identifying corresponding vertices on their convex hull (this leads to four pairs of parallel edges; 
AP^ is pictured in Figure 6.1). The odd Aztec pillowcase OPn is defined analogously using two 
copies of ODn (Figure 6.2 shows OP4; the vertex identifications have been carried out). Note 
that OPn is precisely the planar dual of APn. 

Let pi''' be the graph obtained from the path P„ by weighting each of its edges by q. 

Let Sn be the directed graph obtained from P^2i by regarding each edge as a pair of anti- 
parallel arcs of weight 2, and including an extra vertex v, an arc of weight 2 directed from v to 
the rightmost vertex of P^-i^ ^^nd a loop at v of weight 4 {S3 is pictured in Figure 6.3). Let 5^ 

be constructed analogously from Pn-i , with the loop at the new vertex having weight —4 (^5 is 
illustrated in Figure 6.4). 

30 



Figure 6.2. The odd Aztec pillowcase OP4. 




Figure 6.3. S5. Figure 6.4. S'^. 



Theorem 6.1. AP^ is similar to the disjoint union of two copies of ADn-i, one copy of S2n-i, 
and one copy of S'2n-i- 

Proof. Apply Lemma 1.1 to obtain that APn is similar to the disjoint union of ADn-i and 
the weighted directed graph ADn obtained from ADn by marking the vertices on its convex hull, 

replacing each edge between two unmarked vertices by a pair of anti-parallel arcs of weight 1, 
and replacing each edge between an unmarked vertex v and a marked vertex w by an arc {v, w) 
weighted 1 and an arc {w, v) weighted 2. 

Define the operators N, S, E and W on the set of vertices of ADn to be the nearest neighbors 
in specified cardinal directions, in analogy to the operators NE, NW, SW and SE from the proof 
of Theorem 3.1. Let V be the set of vertices of ADn, and Vi the subset of vertices that are not 
on the boundary of its convex hull. Clearly, the subgraph of ADn induced by Vi of that copy is 

isomorphic to ADn-i- 

For V €Vi, define 

fv ■= C?N(,;)eN(j;) " d'w(v)ew(v) + <^S(?;)es(i;) - C^E(j;)eE(j;) , (6.1) 

where d^ = ^ t£ w G Vi, dw = 2 if w G V \ Vi, and doc = 0. 

A case analysis similar to that in the proof of Lemma 2.3 checks that {fv'vGVi} spans an 
F-invariant subspacc V of U = U^^^ , and that the matrix of the restriction of to V in the 
basis defined by (6.1) is just the adjacency matrix of ADn-i- 

Let Ui , . . . , U2n be the vectors in U that are linear combinations with coefficients or 1 of the 
basis vectors {e„ : v G V}, and whose supports are given by the patterns in Figure 6.5. Let be 
obtained from Uk by multiplying the coefficient of e„ by Cy, where Cy equals 1 for black vertices 

and —1 for whites. 

One readily checks that {ui, . . . ,U2n} are linearly independent; denote their span by Ui. It 
is easy to verify that Ui is F-invariant, and that the matrix of the restriction of F to Ui in the 
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Figure 6.5. The vectors Uk for n = 3. 



basis {ui, . . . ,U2n} is 



2 

2 










2 

2 

2 









(6.2) 



A similar statement holds for the vectors {u'l, 



i2„}, with the restriction of F to the sub- 



space they span being given by the negative of the above matrix. An approach similar to the 
one detailed for three previous examples in this paper shows that the vectors {/„ : v £ Vi}, 
{ui, . . . , U2n} and {u[, . . . , ui^^} span the entire space U, and thus form a basis of it. The matrix 

of F in this basis has the desired block diagonal form. □ 



Corollary 6.2. The number of spanning trees of the Aztec pillowcase graphs is given by 

2 2n-l / j_\ 2 

n 



t(AP„) = t(OP„ 



n 

l<j,k<2n-l 



jir kir \ 
4 cos — cos — . 
2n 



2n J 



(6.3) 



Proof. The first equality follows from the fact that 0P„ is the planar dual of AP„. The spec- 
trum of the matrix (6.2) is readily determined. It follows from Theorem 6.1 and [6, Equations (5)- 
(6)] that 



P(^P„;a;) = (a;-4)(a; + 4) JJ (a; -4 cos 



2n 



n 

l<j,k<2n-l 



ITT klT 

4 cos — cos — 
2n 2n 



(6.4) 



Since all vertex degrees in j4P„ are equal to 4, the eigenvalues of its negative Laplacian are 
obtained by adding 4 to the eigenvalues of the adjacency matrix of ADn, which in turn are im- 
mediately read off from (6.4). Then (6.3) follows by a well known result (see e.g. [7, Theorem VI]) 
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which states that the number of spanning trc;es of a connected graph on m vertices is equal to 
1/m times the product of the non-zero eigenvalues of its negative Laplacian. □ 



Remark 6.3. Embed APn on the pillowcase manifold (topologically equivalent to S*^). Let 
AOPn be the graph obtained by superimposing upon APn its planar dual (which as we pointed 
out is isomorphic to OPn), regarding all points where edges cross or meet as vertices, and all 
lines connecting such points as edges. Then Theorem 4.2 implies that 

M(AOP„ \ {v, w}) = t(AP„) (6.5) 

for any vertex v of APn and any vertex w of OPn on the same face of AOPn as v. By (6.3) we 
obtain a product formula for M(AOP„ \ {u, w}). In particular, this number is independent of v 
and w, a fact not apparent a priori. 

7. Some open problems 

The number of perfect matchings of the (2n-|- 1) x (2n + 1) grid i?„ with a central unit hole does 
not seem to factor into a product of small primes. We have the following prime factorizations: 



M{Hi) 


= 2 




M(i?2) 


= 22. 




M(iJ3) 


= 23. 




M{H4) 


= 2^. 


61212 




= 2^. 


3l2 • 1132 • 27l2 


M{He) 


= 26. 


5924421592 


M{Hr) 


= 2^. 


74172 • 1326051292 


M{Hs) 


= 28. 


44812 • 85132 • 99292 • 163612 


M{Hq) 


= 29. 


46392 • 233576763339021112 


M(Hio) 


= 2^° 


. 72 . 732 . 19l2 . 4792 . 511512 . 29056107452232 


M(ifn) 


= 2" 


• 10332 . 10492 . 16632 • 1661512 • 42412867396854492 


M{Hi2) 


= 2^2 


• 4l2 . 1372 • 70572 . 20992575527970355281835400921 



As pointed out in [6], the factors in a product formula of the kind presented in the previous 
sections can be grouped into small groups so that the product over each group is an integer. This 
shows that such products have a fairly large number of prime factors. In this light, the above 
data suggests that it is unlikely for M(if„) to possess a similar product expression. It would be 
interesting to see whether M(7J„) can nevertheless be expressed in terms of such products, e.g. 
as a sum of a small number of them. 

We conclude with an open problem that suggests a possible generalization of the linear squar- 
ishness of the square grid to grids in higher dimensions. 
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Let Gn ' be the subgraph of the d-dimensional cubical grid induced by the lattice points 
with positive coordinates less or equal to d (thus G„ is just the square grid of Section 2). 
Calculations show that for d = 3 their characteristic polynomials start out as 

P(G'p);a;) =a; 

ViG'i^-x) = (a;-3)(a; + 3) [{x-l){x + l)f 
V{Gf;x) = xix^ - 18) [x^{x^ - 8)(x2 - 2)2]3 
P{Gf^;x) = {x^ - 3a; - 9){x'^ + 3a; - 9) 

X [(a;^ - 3a; + l){x^ +3x + l){x'^ -x- ll)(a;^ + x- ll)(a;2 -x- lf{x'^ +x- iff 
P(Gf ^a;) = xix - 3)(a; + 3){x'^ - 27) 

X [{x - 2){x + 2){x^ - 12)(a;2 - ix + l){x^ + ix + l){x^ - 2x - 11) 

X (a;2 + 2a; - ll)(a;2 - 2a; - 2f{x^ + 2x- 2f{x - lf{x + lf{x^ - Sff. 

This suggests that Gn niay be similar to the disjoint union of three copies of some graph on 
{n^ — n)/3 vertices with some other graph on n vertices. It would be interesting if one could 
prove this, and identify the smaller graphs. 

The pattern for d = 4 is equally suggestive. We have 

PiG^^^;x)=x 

PiG^2^;x) =x'^{x- A){x + 4) [x{x - 2)(x + 2)]^ 

P(G^^); a;) = x'^{x'' - 32)(a;2 - 8)^ [a;*(a;2 - 18)(a;2 - ^fix"" - 2ff 

P(Gf a;) = x^{x - 2f{x + 2f{x'' - 4x - 16){x^ + Ax - 16){x^ - 20f 

X [a;^(a; - 2)(a; + 2){x^ - 20){x^ - 4a; - l){x^ + 4x - l){x^ - 2x - 19){x^ +2x- 19) 
X (a;2 - 2a; - 4)^(a;2 + 2x- 4f{x - lf{x + ifix"^ - 5ff 
P(G^^^ a;) = x^ix - 4)(x + 4)(a;2 - 48)(a; - 2f{x + 2f{x'^ - 12f{x^ -Ax- 8f{x^ + Ax- Sf 



X [x^^{x - 


-3)(a;- 


f 3)(a;2 - 6a; + 6){x'^ + 6x + 6)ix^ - 2x - 


-26)(a;2 + 2a;-26)(a;2-27) 


X 


- 4a; - 


8)(a;2 + 4a; - 8){x'^ - Ax + ifix"^ + Ax - 


H)3 


x(a;2 


-2a;- 


ll)=^(a;2 + 2a; - ll)^(a; - 2f{x + 2f{x^ 


- 12)5 


X {x^ 


-2x- 


2)9(a;2 + 2x- 2f{x - lf\x + l)i°(x2 


-3fX. 



This raises the problem of understanding whether G„ is similar to the disjoint union of four 
copies of some graph on (n^ — n^)/4 vertices with some other graph on vertices. More generally, 
it would be interesting to see if an analogous phenomenon occurs for general d. 

Another open problem (suggested by one of the referees) is to consider mixed Aztec pillowcase 
graphs and see if they admit a decomposition similar to the one given in Theorem 6.1. 
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